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Abstract 

A kinetic approach to the Bose-Einstein condensates (BECs) is proposed, based on the Wigner- 
Moyal equation (WME). In the quasi-classical limit, the WME reduces to the particle number 
conservation equation. Two examples of applications are: i) a self-phase modulation of a BE 
condensate beam where we show that a part of the beam is decelerated and eventually stops as a 
result of the gradient of the effective self-potential; ii) the derivation of a kinetic dispersion relation 
for sound waves in the BECs, including a collisionless Landau damping. 
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I. INTRODUCTION 



Presently, the Bose-Einstein condensates (BECs) provide one of the most active and 
creative areas of research in physics [1, 2]. The dynamics of the BECs are usually described 
by a nonlinear Schrodinger equation (known in this field as the Gross-Pitaevskii equation 
(GPE) [3, 4]), which determines the evolution of a collective wave function of ultra-cold 
atoms in the BECs, evolving in the mean field self-potential. 

In this paper, we propose the use of an alternative but nearly equivalent approach to the 
physics of BECs, based on a kinetic equation for the condensate. We also show that this 
kinetic theory can lead to a more complete understanding of the physical processes occurring 
in the BECs, not only by providing an alternative method for describing the system, but 
also by improving our global view of the physical phenomena. It is our hope that this will 
also lead to the discovery of new aspects of BECs. 

The key point of our present approach is the use of a Wigner Moyal equation (WME) 
for the BECs, describing the spatio-temporal evolution of the appropriate Wigner function 
[5]. Wigner functions for the BECs were discussed in the past [6, 7] and the WME has been 
sporadically used [8]. But no systematic application of the WME to BECs has previously 
been considered. In the quasi-classical limit, this equation reduces to the particle number 
conservation equation, which is a kinetic equation formally analogous to the Liouville equa- 
tion, but with a nonlinear potential. A description of the BECs in terms of the kinetic 
equation is adequate to deal with a series of problems, as exemplified here, and can be seen 
as an intermediate (in accuracy) between the GPE and the hydrodynamic equations usually 
found in the literature. 

The manuscript is organized in the following fashion. In Section 2, we establish the WME 
and discuss its approximate version as a kinetic equation for the Wigner function. We then 
apply the kinetic equation to two distinct physical problems. The first one, considered in 
Section 3, is the self-phase modulation of a BEC beam. A similar problem has been studied 
numerically in the past [9]. Here, we derive explicit analytical results and show that a part 
of the BEC beam is decelerated and eventually comes to a complete halt as a result of the 
collective forces acting on the condensate. The second example is considered in Section 4, 
where we establish a kinetic dispersion relation for sound waves in the BECs, giving a kinetic 
correction to the usual Bogoliubov sound speed [12, 13] and predicting the occurrence of 
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Landau damping [10, 11]. Our description of Landau damping is significantly different from 
that previously considered for transverse oscillations of BECs [14]. Finally, in Section 5, the 
virtues and limitations of the present kinetic approach are briefly discussed. 



II. WIGNER MOYAL EQUATION FOR THE BOSE CONDENSATE 

It is known that for an ultra-cold atomic ensemble, and in particular for BECs, the ground 
state atomic quantum field can be replaced by a macroscopic atomic wave function ip. In a 
large variety of situations, the evolution of ip is determined by the GPE 

^ = -^ V V + (^o + V e// M (1) 

where Vo = Vo(f) is the confining potential and V e ff is the effective potential which takes into 
account the inter-atomic interactions inside the condensate, as determined in its simplest 
form by V e ff(f,t) = g\ip(r,t)\ 2 , where g is a constant [3, 4]. 

Let us consider the situation where this wave equation can be replaced by a kinetic equa- 
tion. In order to construct such an equation, we introduce the Wigner function associated 
with ip, such that [5] 

W(r, k,t) = J ip(r + s/2, t)ijj*(f- s/2, t) exp(-k ■ s) ds. (2) 

It is then possible to derive (see the Appendix) the following evolution equation for the 
Wigner function 

(—k -V -ih^-]w = -2VsinAW, (3) 
\2m dtj w 

where a bi-directional differential operator is denoted by 

This operator acts to the left on V and to the right on W [5]. In this equation, the potential 
is 

/— * d k 
W(r,k,t)— + 5V, (5) 

where 
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5V = g^t)\ 2 -jw^k,t)j^j (6) 

can be considered a noise term associated with the square mean deviations of the quasi prob- 
ability, determined by the Wigner function W with respect to the local quantum probability, 
determined by the wave function ip. 

Equation (3) can be seen as the WME describing the space and time evolution of the 
BECs, and it is exactly equivalent to the GPE (1). However, it is of little use in the above 
exact form, and it is convenient to introduce some simplifying assumptions. This is justified 
for the important case of slowly varying potentials. In this case, we can neglect the higher 
order space derivatives, and introduce the approximation sin A ~ A. This corresponds to 
the quasi-classical approximation, where the quantum potential fluctuations can also be 
neglected, viz. 5V — > 0. Introducing these two simplifying assumptions, valid in the quasi- 
classical limit, we reduce the WME to a much simpler form 

where v = hk/m is the velocity of the condensate atoms corresponding to the wavevector 
state k, and F = W is a force associated with the inhomogeneity of the condensate self- 
potential. The nonlinear term in the GPE (1) is hidden inside this force F. As we will 
see, this nonlinear term will look very much like a ponderomotive force term, similar to the 
radiation pressure. 

It should be noticed that this new equation is a closed kinetic equation for the Wigner 
function W . In this quasi-classical limit, W is just the particle occupation number for 
translational states with momentum p = fik. Equation (7) is equivalent to a conservation 
equation, stating the conservation of the quasi-probability W in the six-dimensional classical 

— * 

phase space (r, k), and can also be written as 

jW(r,k,t) = 0. (8) 

This kinetic equation can then be used to describe physical processes occurring in a BEC, 
as long as the quasi-classical approximation of slowly varying potentials is justified. The 
interest of such a kinetic descriptions will be illustrated with the aid of two simple and 
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different examples, to be presented in the next two sections. Many other applications can 
be envisaged, and will be explored in future publications. 



III. SELF-PHASE MODULATION OF A BEAM CONDENSATE 

Let us first consider the kinetic description of self-phase modulation of a BEC gas, moving 
with respect to the confining potential Vo(f). Here, we can explore the analogy of this 
problem with that of self-phase modulation of short laser pulses moving in a nonlinear 
optical medium and well known in the literature [15]. In order to simplify our description, 
we consider the one-dimensional problem of a beam moving along the z-axis and neglect the 
axial variation of the background potential dVo/dz ~ 0. The radial structure of the beam 
can easily be introduced later, and will not essentially modify the results obtained here. The 
kinetic equation (7) can then be written as 

with v z and F z given by, respectively, 

hk d w „ dk d w , . 

Vz = — + 9 7 rl(z,t , F z = — = -g—I(z,t), 10 
m at dt oz 

where we have used the intensity of the beam condensate, as defined by 

I(z,t) = J W(z,k,t)^. (11) 

Let us assume that the ultra-cold atomic beam has a mean velocity v = hko/m. This 
suggests the use of a new space coordinate rj — z — v t. In terms of this new coordinate, the 
quasi-classical equations of motion of a cold atom in the beam can be written as 



dri dh 1 , s „. 

dk dh g d , 

dt = ~fr, = ~hd~r ] iV,th 

where we have introduced the Hamiltonian function 

k ( k \ g 

Kv, k, t) = u(rj, k, t) - kv = — - - k + -I(r], t). (13) 

m \ 2 / n 
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Here cu(r),k,t) is the Hamiltonian in the rest frame expressed in the new coordinate. A 
straightforward integration of the equations of motion leads to 



k{t)=k,- 9 -f^I{r 1 ,t')dt'. (14) 

At this point it is useful to introduce the concept of the beam energy chirp, < e(r),t) >, 
in analogy with the frequency chirp of short laser pulses [15]. By definition, it will be the 
beam mean energy at a given position and at a given time 

/dk 
W( V ,k,t)u( V ,k,t)—, (15) 

where the weighting function W(rj, k,t) is the solution of the kinetic one-dimensional equa- 
tion (9). A formal solution of this equation can be written as 

W( V , k, t) = W{ Vo { V , k, t), k { V , k, t),t ), (16) 

where r] and k are the initial conditions corresponding to the observed values at time t, as 
determined by the dynamical equations (13). Replacing it in equation (14), we obtain 



< e(rj : t) >= h J W(r] ,k ,t ) 



k 2 g w % 



— . (17) 



From equation (14) we notice that dk = dko- Neglecting higher order nonlinearities, we 
can then rewrite the above expression as 

< e( V , t) >=< e(0) >-^gf £-I( V , t')dt', (18) 

m Jo or] 

where < e(0) >=< e(i] ,to) > is the initial beam energy chirp. 

Let us first consider that the beam profile I{rj) is independent of time. This is, of course, 
only valid for very short time intervals where the beam velocity dispersion is negligible. In 
this simple case, we have 

< e(r], t) >=< e(0) > -%v Q g (^j t. (19) 

The maximum energy shift will be attained at some position inside the beam profile, 
V = Vmax, determined by the stationarity condition 
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d . . (d 2 i\ 

^ <e( ^ } >= [w) = 



(20) 



In order to deduce more specific answers, let us assume a Gaussian beam profile 



/(??) = J exp (W> 2 ), 



(21) 



where a determines the beam width. For this profile we have r\, 
to the following value of the maximum energy shift 



'max 



±a/y/2, which leads 



Ae(i) =< e(t) > 



max 



< e(0) >= ± gv I e 



a 



(22) 



This is similar to a well known result in nonlinear optics, stating that the maximum energy 
chirp due to a self-phase modulation is proportional to t, or to the distance travelled by the 
beam, d = v t. This result clearly indicates that the initial beam will eventually split into two 
parts, one being accelerated to higher translational speeds, and the other being decelerated. 
This would correspond to the red-shift and to the blue-shift observed in nonlinear optics. 
The decelerated beam will eventually stop after a time t ~ r, such that Ae(r) =< e(0) >. 
This will determine the condition for translational beam freezing. 

It should be noticed that the same result could also be obtained directly from the GPE 
(1). But the interest of the present derivation is that it demonstrates the irrelevance of the 
phase of the wave function ip, because such a phase was ignored in our kinetic calculation. 
Therefore, instead of the self-phase modulation, it would be more appropriate to call it a 
beam self-deceleration. 

Another interesting aspect of our kinetic approach is that it can be easily refined, as 
shown briefly here. Let us then improve the above calculation by considering the beam 
dispersion. This will inevitably occur because of the linear velocity dispersion of the atomic 
beam. Such a dispersion will decrease the chirping effect, because of the decrease in time 
of (dl/drj). In order to model it, we can assume a time- varying Gaussian beam shape, as 
described by 



If we now consider that a(t) = cr (l + 5t 2 ), where 5 = (2m/h 2 )Ae /a is proportional to the 
initial energy spread Aeo- This will lead to a maximum energy shift 
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(23) 



Ae d {t) = -^Ae(f). (24) 

It is clear that the linear beam velocity dispersion will decrease the maximum attainable 
chirp, by changing the linearity with time into a logarithmic law. However, this will only 
occur for very long times, t ~ l/VS, which for ultra-cold atomic beams with a very low 
translational energy dispersion Aeo will not be relevant. 

The other cause of the beam dispersion is the nonlinear process itself, which will eventu- 
ally break up the initial pulse into two distinct pulses. In this case, the self-phase modulation 
process will not be attenuated because the beam width will be conserved, but the two sec- 
ondary pulses will suffer self-phase modulation themselves, and will eventually break up 
later, resulting in the formation of several secondary pulses with different mean energy. 
However, the nonlinear dispersion will also be negligible as far as a 2 > Am\Ae(t)\t 2 /h 2 . 
A more complete description of all these dispersion regimes can be obtained by solving 
numerically the kinetic equation (9). 



IV. KINETIC DESCRIPTION OF BOGOLIUBOV OSCILLATIONS 

The second example of an application of the kinetic equation for the BECs deals with the 
dispersion relation of sound waves. For simplicity, we consider again the one-dimensional 
model and neglect the radial structure of the oscillations. This allows us to treat the lowest 
order oscillating modes of the condensate. We assume some given equilibrium distribu- 
tion Wq(z, k,t), corresponding, for instance, to the Thomas-Fermi equilibrium solution in a 
given confining potential Vo(f±, z) [16], and after linearization of the one-dimensional kinetic 
equation (9) with respect to the perturbation W, we obtain 

dt + Vz d~z) W{Z > k > f) + P dk Wo{Z > k > t] = °' (25) 
where the perturbed force is determined by 

Let us now assume perturbations of the form W,I ~ exp(ikz — iut). From the above 
equations we then obtain a relation between the perturbation amplitude of the Wigner 



8 



function W and the perturbed beam intensity / 

where we now specify the particle wavenumber state with k', in order to avoid confusion with 
the wavenumber k of the oscillation that we intend to study. The velocity corresponding to 
this particle state is v' = hk'/m. Integration over the momentum spectrum of the particle 
condensate will then lead to the following expression 

9 r dWo(k>)/dV dk'_ 
1 + h k J (u — hkk'/m) 2ir ~ { ^ } 

This is the kinetic dispersion relation for axial perturbations in the BECs. Let us illustrate 
this result by considering a simple case for a condensate beam with no translational disper- 
sion, or with a translational temperature exactly equal to zero. The equilibrium state of the 
beam can then be described by 

W (k') = 2-irn 5(k' - k' ), (29) 

where no — J Wo(k')dk' /2ir is the particle number density in the condensate. Replacing this 
in the dispersion relation (28). we have 

1 - ~7TTT2 = 0. (30) 

where v' = hk' /m = p' /m is the beam velocity. This can also be written as 

{u - kv' ) 2 = k 2 c 2 s , (31) 

where 



y/gno/m (32) 



is nothing but the Bogoliubov sound speed. Obviously, equation (31) is the Doppler shifted 
dispersion relation of sound waves in the BEC gas, In its reference frame it reduces to 
uj = kc s . 

Let us now consider a situation where, instead of the distribution (29), we have a beam 
with a small translational velocity spread, such that the number of particles with a velocity 
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v' ~ c s is small but nonzero. In this case, the resonant contribution in the integral of equation 
(28) has to be retained, although it is still possible to neglect the kinetic corrections in the 
principal part of the integral. The dispersion relation can then be written, in the condensate 
frame of reference, as 



where k' s = mc s /h is the resonant momentum. The imaginary term in this equation can 
then lead to damping of the sound waves. Writing now uj = kc s + ry, with |7| <C kc s , we 
can then obtain the expression for the damping coefficient 



The above expression corresponds to the non-collisional Landau damping of Bogoliubov 
oscillations in the BECs. The present approach can also be generalized in a straightforward 
way to higher order oscillations of the condensate, where the radial structure has to be taken 
into account [13, 17]. 

V. CONCLUSIONS 

In this paper, we have proposed a kinetic view of the Bose-Einstein condensate physics, 
based on the Wigner-Moyal equation. In the quasi-classical limit, the latter can be reduced 
to a closed kinetic equation for the corresponding Wigner function. The kinetic approach 
to BEC can be seen as an intermediate step between the GPE and the hydrodynamical 
equations for the condensate gas, often found in the literature. 

We have discussed two different physical problems, in order to illustrate the versatility 
of the kinetic theory. One is a self-phase modulation of a BEC beam. The other is the 
dispersion relation of the Bogoliubov oscillations in the condensate gas. The first example 
shows that, due to the influence of its own inhomogeneous self-potential, nearly half of 
the beam is accelerated while the other half is decelerated. Under certain conditions, the 
decelerated part of the beam will tend to a state of complete halt. The second example 
shows that a kinetic dispersion relation for sound waves in BECs can be established, where 
Landau damping is automatically included. The present results only considered the lowest 




(33) 




(34) 
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order modes, but the same approach can be used to describe higher order oscillations of the 
BECs, including their radial structures, as well as the coupling with a background thermal 
gas. This investigation is, however, beyond the scope of the present work. 

Several other different problems relevant to BECs can also be considered in the frame of 
the kinetic theory, such as modulational instabilities [18] and the wakefield generation. This 
indicates that the kinetic theory is a very promising approach to the physics of BECs, which 
will eventually allow us to introduce new ideas in this stimulating area of research, and to 
suggest new configurations to the experimentalists. However, the present work also clearly 
states that the present theory is only valid in the quasi-classical limit and for that reason 
some relevant problems, where the phase of the BEC wave function plays an important role, 
can only be treated by means of the GPE. Surprisingly, the self-phase modulation is not one 
of them, as demonstrated here. 
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APPENDIX: Derivation of the Wigner-Moyal equation 

In the present derivation we follow a procedure already used in other cases, for instance, in 
the case of electromagnetic waves moving in a space and time dependent dielectric medium 
[19]. For a different but nearly equivalent derivation of the WME see the appendix of 
reference [8]. Let us consider two distinct sets of values for space and time coordinates, 
(r*i, ti) and (r 2 , t 2 ), and let us use the notation ipj = ip(fj, tj) and Vj = V(fj, tj), for j = 1,2. 
This allows us to write two versions of the GPE (1) as 

Multiplying the equation j = 1 by ip 2 and the conjugate of the equation j — 2 by tpi, and 
subtract the resulting equations, we obtain 



2m y 1 21 \dt 1 dt 2 



C12 = -(^i - V 2 )C 12 , (36) 



where we have used C\ 2 = ipiip 2 . The above equation suggests the use of two pairs of space 
and time variables, such that 



f 1= f-s/2 , t 1 = t-r/2, 
r 2 = r + s/2 , t 2 = t + t/2. 



(37) 



We can then rewrite the above equation as 



v h 2 d d .d 
- — in- 



C 12 = -(V 1 -V 2 )C 



12- 



m df ds dt 

It can also easily be shown, by developing the potentials Vj around V(f,t), that 



(y 1 -y 2 ) = 2sinh(f.^ + l|)n^). 



(38) 



(39) 



Let us now introduce the double Fourier transformation of the function C± 2 = C(f, s, t, r) 
on the variables s and r, as defined by 



W(f, t, u, k) — J ds J drC(f, s,t,r) exp(—ik • s + iur). 
This can be rewritten in terms of the wave function ip as 



(40) 
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W(f, t, u,k) = J dsj dr^(f + s/2, t + r/2)^*(f - s/2, i - r/2) exp(-i£; • s + iur). (41) 
Using such a development in equation (38), we obtain for the Fourier amplitudes 

[ —k ■ — - h— W — —2V sin A'W, (42) 
\m or at J 

where we have used the following differential operator, operating backwards on the potential 
V(r,t) and forward on W, viz. 

A '=2 (5f- ag- Sac] ' (43) 

This is a formidable equation for W, which can be simplified by noting that the GPE 
implies the existence of a well defined relation between energy and momentum. This means 
that uj must be equal to some function of k, or u = cu(k). Hence, we can state that 

W(f, t, uj, k) = W(f, k, t)5(co - u(k)). (44) 

This leads to a much simpler evolution expression for W(f, k, t). Before writing it down, 
we should also notice that the nonlinear term in V depends on \ip\ 2 , and not on the function 
W. Thus, we can finally write 

' V " iH o\) W = ~ 2{Vo + 9\^\ 2 )(^ A W, (45) 
where A is a simpler differential operator defined by 

— * 

The function W(f, k, t) can be seen as the Wigner function associated with the GPE, 
and equation (45) as the WME equation that describes its spatio-temporal behavior. This 
equation is equivalent to the initial wave equation (1), but it is not a closed equation for 
the quasi-probability function W. Therefore, some simplifying assumptions have to be 
introduced in order to make it more tractable, as explained in Section 2. 
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